The BLAST Survey of the Vela Molecular Cloud: 
Dynamical Properties of the Dense Cores in Vela-D 
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ABSTRACT 

The Vela-D region, according to tlie nomenclature given by iMurphv fc Mavl (|l991l ). of the 
star forming comp l ex kn own as the Vela Molecular Ridge (VMR), has been recently analyzed 
in details by lOlmil ( 20091 ) . who studied the physical properties of 141 pre- and proto-stellar cold 
dust cores, detected by the "Balloon-borne Large-Aperture Submillimeter Telescope" (BLAST) 
during a much larger (55 deg^) Galactic Plane suryey encompassing the whole VMR. This survey's 
primary goal was to identify the coldest, dense dust cores possibly associated with the earliest 
phases of star formation. In this work, the dynamical state of the Vela-D cores is analyzed. 
Compari son to dyna mical masses of a sub-sample of the Vela-D cores estimated from the ^^CO 
survey of lElial ( 20071 ), is complicated by the fact that the -'^'^CO linewidths are likely to trace the 
lower density intercore material, in addition to the dense gas associated with the compact cores 
observed by BLAST. In fact, the total internal pressure of these cores, if estimated using the 
^■^CO linewidths, appears to be higher than the cloud ambient pressure. If this were the case, 
then self-gravity and surface pressure would be insufhcient to bind these cores and an additional 
source of external confinement (e.g., magnetic field pressure) would be required. However, if one 
attempts to scale down the "'^^CO linewidths, according to the observations of high-density tracers 
in a small sample of sources, then most proto-stellar cores would result effectively gravitationally 
bound. 

Subject headings: submillimeter: ISM — stars: formation — ISM: clouds — balloons 
1. INTRODUCTION 
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Pre-stellar cores represent a very early stage 
of the star formation (SF) process, before col- 
lapse results in the formation of a central pro- 
tostar. The physical and dynamical properties 
of these cores can reveal important clues about 
their nature. More than a thousand new pre- 
and proto-stellar cores have been identified by 
BLAST, the Balloon- borne Large-A perture Sub- 
millimeter Telescope ( Pascald 120081 ) . during its 
second long duration balloon (LDB) science flight 
in 2006. BLAST detected these cold cores in 
a '^ 50 dcg^ map of th e Vela Molecular Ridge 
(VMR) (Ncttcrfieldll2009l ). and was capable of con- 
straining the temperatures of the detected sources 
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using its three-band photometry (250, 350, and 
500 /im) near the peak of the cold core SED. 



The VMR (jMurphv fc Mavlll99lL iLiseau et al. 



19921 ) is a giant molecular cloud complex within 
the galactic plane, in the area 260° < / < 272° 
and —2° ^6^3°, hence located outside the 
solar circle. Its m ain propert i es hav e be en re- 
cently revisited by iNetter field (2009[) and lOlmi 



(2009 1. In particular, [Ohni (2009i ) used both 



BLAST and archival data to determine the spec- 
tral energy densities (SEDs) and the physical pa- 
rameters of each source detected by BLAST in 
the smal ler region of V ela-D, where observations 
from IR (lGia.nnini 20071 ) to millimeter wavelengths 
( Massi et al.ll2007 ) were already available. 

In this wor k we will use the spectral line data 
of lElial ([20071) to perform an analysis of the dy- 
namical state of the cores. In Section [2] we discuss 
which cores are gravitationally bound. Then in 
Section |3] we discuss the effects of pressure. We 
draw our conclusions in Section |4l 

2. DYNAMICAL STATE OF THE CORES 
2.1. Core Size Distribution 

The cores in Vela-D have a similar range of sizes 
((i?^) = 0.09±0.02 pc) to those in the overall VMR 
{{Re) = 0.07 ± 0.02 pc), for comparable range of 
masses. Please note that although we use the term 
core size, or core radius, what we actually measure 
with BLAST is an "effective core width", set by 
the power-law distribution of the column density 
in the core and by the BLAST beam. The small 
range of sizes for a large range of masses in both 
VMR and Vela-D is similar to what Enoch et al. 
(2008) found for the proto-stellar cores in Serpens, 
Ophiucus and, to a lesser degree, in Perseus. 

We also evaluated the ratios of angular de- 
convolved size (from BLAST) to beam size, 
^doc/^mb- Power-law density profiles in pre- 
stellar cores have received both a theoretical 



(see for example the discussions in lYoung et al 
20031 and McKee fc Ostriker 2007) and observa- 
tional support (e. g., IWard-Thompson et al.lll999l 
Andre et al.ll20Q0l and references therein). Assum- 
ing that also the Vela-D cores can be described 
by power- law dens ity profiles, then according to 
Young et al.l (|2003l ) ^doc/^mb should depend only 
on the index of the power law, and not on the 
distance of the source. We find ^dcc/^mb to have 
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Fig. 1. — Gray-scale image shows the BLAST 
250 /xm map of Vela-D, with galactic coordinates 
in degrees. Solid c ontours re present the ^■^C0(2 — 
1) emission (from lElial 120071 ). integrated between 
2 km s~^ and 15 km s~^. The first contour is at 
2Kkm s^^, and contour spacing is ~ 7Kkm s^-'^. 



a median value of 1.41 for Vela-D, correspond- 
ing to a density power-law exponent p ~ 1.6 
(JYoung et al.ll2003l ). This value for t he power-law 
expon ent, similar to that found by Enoch et al. 



(|2008[ ) in Ophiucus, is consistent with mean p- 
values found from radiative t ransfer modeUng o f 
Class and Class I envelopes (JYoung et al.ll2003l ) . 



2.2. Bonnor-Ebert Analysis 

In this section we investigate the dynamical 
state of starless and proto-stellar cores, in order 
to determine whether or not the BLAST starless 
cores in Vela-D are pre-stellar, i.e. they will col- 
lapse and form one or more stars. We first analyze 
the dynamical state of the cores by calculating the 
critical Bonnor-Ebert mass, Mbe- In this analysis 
we assume that all cores are spherical in shape, as 
our angular resolution at the distance of Vela-D 
is not sufficient to determine with enough accu- 
racy the ellipticity of the cores. Because elongated 
cores can be considerably more stable than ex- 
pecte d from a pure spheric al Bonnor-Ebert model 
(e.g.. lStahler fc Pallall20041 ) . then we may be over- 
estimating the non-stability of some cores. Fur- 
thermore, it should be noted that even if the core 
mass is larger than the Bonnor-Ebert mass, this 
condition alone does not guarantee that the core 
is truly gravitationally unstable, as we are neglect- 
ing any other non-thermal support. Other sources 
of support are discussed in the next sections. 

The Bonnor-Ebert mass, Mbe, can be deter- 
mined from the observables, i.e. the core radius 
and temperature: 



MBE[Mo] = 1.8i?BEMT[K] 



(1) 



where T is the core temperature, as determined 
from the SED fits (,OlmL2009i ). and i?BE represents 
the core radius at which the internal and ambient 
pressures are equal, for a Bonnor-Ebert sphere in 
hydrostatic equilibrium. 

i?BE can be estimated as Rbe = f Rdcc, where 
i?dcc is the measured linear deconvolved half-width 
at half-maximum and / is a numerical factor that 
has been determined as follows. For each core 
the approximated Bonnor-Ebert density profile for 
an isothermal sphere, flat at the core centeilll and 
falling off with r^^ to the outside (see Figure [2]), 
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Fig. 2. — Example showing how i?BE is esti- 
mated. The red solid line represents the profile 
of the BLAST 250 /im beam at the distance of 
Vela-D. The blue solid line shows the approxi- 
mated Bonnor-Ebert density profile, from which 
the simulated i?BE is calculated. The black solid 
line represents the convolution of the beam and 
density profiles, and the dashed black line shows 
the Gaussian best-fit to the convolved profile, from 
which the source simulated radius, R^^"^ (HWHM), 
is determined. 



^The core center is defined here as the region where the 
normaUzed radius r'n < 1, with r'n = r <^ iTTGrn^cv po / (kT) , 



is generated. The density profile is obtained using 
the temperature and average density estimated for 
each core, and the theoretical radius, Rbe, of this 
sphere is estimated. This profile is then convolved 
with the BLAST beam at 250 /im and then fitted 
with a Gaussian, obtaining a simulated "observed" 
radius i?*''™ (defined as the half-width at half- 
maximum of the fitting Gaussian). Finally, this 
simulated radius is deconvolved, obtaining i?!™, 
and then compared with Rbe- Their ratio is then 
the factor / = -RBE/^dcc specific to that source. 
The factor / can then be used to determine Rbe 
from the measured radius, i?dec, thus keeping into 
account any mismatch between i?|™ and i?doc. 

At the distance of Vela-D, and for the range of 
core temperatures found in this cloud, we find an 
average value (/) = 2.55 ± 0.14. Thus, from the 
observed i?doc and core temperature, the value of 
Mbe can be estimated. To further rid the sam- 
ple of unrealistically small core sizes we also elimi- 
nated all data for which the beam rad ius exceeded 
the deconvolved source radius. See Ijiiina et al. 



(|l999l ) and lEnochet all (|2007[ ) for further discus- 
sion on the effects of instrumental resolution. 

Figure [3] plots the total core mass, Mcorc, also 
determined from the SED fits, vs. the Bonnor- 
Ebert mass for all BLAST cores in Vela-D that 
have well-determined SED fits. We also plot the 
completeness limits and a vertical -dashed line cor- 
responding to the estimated minimum measurable 
Bonnor-Ebert mass, calculated as; 

M^i^[MQ] = l.SR^tlvc] Tcompi[K] (2) 

where Teompi = 12 K, and i?ff = (/}i?bcam, 
where i?bcam corresponds to the beam linear half- 
width at half-maximum at the distance of Vela-D. 
The cores with Mcore > Mbe represent sources 
which, in the absence of further internal support 
(turbulent, magnetic, etc.), are gravitationally un- 
stable and will collapse (or are already collapsing) 
in the presence of external perturbations. With- 
out further knowledge of the internal structure of 
these cores we cannot yet establish whether they 
are collapsing or will collapse to form stars. On 
the other hand, the cores with Mcorc < Mbe rep- 
resent sources which, assuming only thermal sup- 
port as in Bonnor-Ebert spheres, are stable against 
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Fig. 3. — Total core mass, derived from the SED 
fits, vs. the Bonnor-Ebert mass. Red symbols are 
for starless cores, black symbols for proto-stellar 
cores. The median value of Mbe is 5.7 M©. The 
solid line corresponds to Afcorc = Mbe- The two 
horizontal lines correspond to the completeness 
limits for sources with T > 12 K (dotted line) 
and T > lOK (dashed line), respectively (jOlmi 
20091 ). The vertical dot-dashed line correspond to 
the estimated minimum measurable Bonnor-Ebert 
mass, M^^^. 



where mav is the mass of the average gas particle and po 
is the density at the center of the core. 



external perturbations. Any further internal sup- 
port will thus add to their stability against gravi- 
tational collapse, and might eventually lead to the 
core disruption. 

We find a median McorcZ-'^^BE = 0.9 for all 
sources in Vela-D, i.e. very near the critical value, 
and indicating that we have about the same num- 
ber of sources with Mcore < Mbe and Mcore > 
Mbe- However, if we restrict ourselves to mass 
ranges for which we are progressivel y more com - 
plete (Mcoro > 4Mq, see Table 3 of lOlmil [2OO9I) . 
our sources exceed their corresponding Bonnor- 
Ebert mass by a median factor of A/core /-^j^be = 
1.5 — 3.3. Figure [3] shows very little difference be- 
tween starless and proto-stellar cores, and we can 
also see that the core mass range where the critical 
value Mcoro = Mbe is crossed is Mcore ~ 3 — 8Mq. 
In this range we are only partial ly complete (see 
NetterfieldlT2009l and lOlmil l2009l ) and thus more 
sensitive measurements are needed in order to con- 
firm that this is an evolutionary effect. 

2.3. Virial Analysis 

2.3.1. Determination 0/ Virial Masses 

We now investigate the dynamical state of the 
Vela-D cores by estimating their virial masses. 
The comparison of dust-derived masses to molecu- 
lar line observations gives a more robust method of 
determining if the cores are gravitationally bound. 
In fact, while our determination of the dust masses 
provide an estimate of the gravitational potential 
energy, the molecular linewidths constitute a mea- 
sure of the internal energy of the cores, provided 
that the selected molecular lines are truly trac- 
ers of the cores observed at (sub)millimeter wave- 
lengths. This is typically the case with molecu- 
lar probes of the interior of pre- and proto-stellar 
cores, such as NH3 and N2H+, whereas more abun- 
dant and diffuse molecules, such as ^'^CO described 
below, will provide only an upper limit to the in- 
ternal energy of the cores. 

Here we wi ll be using the ^^C0(2 — 1) data 
of lElial (|2007l ). but instead of using the cata- 
log of these au t hors (obtained by CLUMPFIND, 



Williams et al.l 119941) to associate a molecular 



clump to a BLAST core, we use the coordinates 
of the BLAS T objects (when they fall within the 
^^CO map of lEliall200lT ) to find the corresponding 
"'^^00(2 — 1) spectrum along the line of sight. This 



procedure has two advantages: first, we avoid as- 
sociating the same catalog clump (generally some- 
what more massive and larger in size) to multiple 
BLAST cores, and in addition we can directly per- 
form a Gaussian fit of the ^^C0(2 — 1) line without 
having to rely on the ability of CLUMPFIND to 
separate adjacent (in space and velocity) clumps. 

We selected only BLAST cores with corre- 
sponding well-determined ^■^C0(2 — 1) linewidths 
(a total of 35 sources). In the spectra of four 
sources we find multiple ^'^CO velocity compo- 
nents along the line-of-sight to the BLAST core. 
However, in three of these sources we observe only 
two velocity components, and their linewidths are 
very similar or they differ by less than 50 %. Thus, 
as long as the linewidths are not too different, it 
does not matter which of the two -'^'^00(2 — 1) ve- 
locity components is actually associated with the 
BLAST core. In the remaining source we find mul- 
tiple velocity components with linewidths ranging 
from 0.3km s~^ to 1.8km s^^, and because we 
cannot know with certainty which one is actually 
associated with the BLAST core, we did not use 
it for the determination of the virial mass. 

We estimate the virial mass of the cores as 



( MacLaren et al.lll988D : 



Mvir[Mo] = 210i?dcc[pc] (Ay[km s-^]f (3) 

where ISV represents the total linewidth of the 
mean particle, and is calculated as the sum of the 
thermal and turbulent components: 



Ay2 = Ay^^gco 



fcr81n2 



1 



1 



(4) 
where IWi^co is the "C0(2 - 1) line FWHM, 
ruav = 2.3 amu is the mean molecular weight and 
7TT,i3CO = 29 amu is the molecular weight of ^^CO. 
Figure m plots the total core mass Mcore vs. the 
virial mass Mvir for all the BLAST cores in Vela- 
D selected as described above. The solid line in 
this figure represents the self-gravitating limit de- 
fined by K = —U, where K is the kinetic and 
U the gravitational potential energy, correspond- 
ing to Mcoro — 0.5Afvir- We note that one could 
expand the definition of Bonnor-Ebert mass to in- 
clude the effects of turbulence, thus writing Eq. ^ 
as: 

Mbet[M0] = 1.8i?BE[pc] Teff [K] (5) 
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Fig. 4. — Total core mass, derived from the SED 
fits, vs. the virial mass, calculated using the veloc- 
ity linewid th derived from the ^"^00(2 — 1) data of 
Elial (J2007() . The solid line indicates the minimum 
■^core = O.SMvir = O.SA/bet for which the cores 
should be self-gravitating. Black asterisks repre- 
sent the results when no scale factor is applied to 
the ^^C0(2 — 1) linewidths, whereas the red open 
squares show the resulting parameters when the 
linewidths are scaled (with the 1.6 factor, see text) 
and thus are thought to represent more closely the 
core internal velocity dispersion. For the sake of 
simplicity, both starless and proto-stellar cores are 
represented by the same symbols. 



where the effective temperature. Toff, is defined 
as the temperature that describes the total line 
width, i.e.: 



AV'^ 



fcr^ff81n2 



(6) 



= 50Ay2[kms-i]. 
Mbet[Mo] = 90/i?dcc[pc] AV^ [km s-^] , (7) 



which also gives Teff [K 
Eq. (O then becomes: 



where we have used i?BE = / ^doc (see Sec- 
tion [521) • We therefore see that A/bet — M^it, 
and the self-gravitating limit corresponds to 

Mcore = 0.5 Mvir ^ 0.5 MbET- 

In Figure H] we can see that nearly all of the 
cores lie below the self-gravitating line, indicating 
that they are unlikely to be gravitationally bound, 
in agreemen t with the findings of Elia (2007 ) and 
Massi et al.l (|2007t ). However, the i3cO(2 - 1) 
emission may be partially optically thick along 
the line o f sight of the ^'^CO clump centroid (see 
Elial 120071 ) ■ and even when it is indeed optically 
thin it traces all molecular gas along the line of 
sight. Although, as we discussed earlier, we mea- 
sure the ^■^CO linewidths along the line of sight 
to the BL AST c ores a nd not toward the clumps 
centroid of lElial ( 20071 ). we are likely to trace the 
lower density intercore material in addition to the 
dense gas associated with the compact cores ob- 
served by BLAST. We attempt to correct for this 
effect in the next section. 

2.3.2. Application of a Scaling Factor for the 
Linewidths 

Recent follow-up observations toward a sam- 
ple of the Vela-D cores (Olmi et al., in prepara- 
tion) have allowed us to determine the average 
of the i3CO(2 - 1) to N2H+(1 - 0) linewidth ra- 
tio, equal to ~ 1.6, in a small sub-sample of 7 
cores. For comparison, in the Pipe cores Rath- 
borne et al. (2008) find that the C^ ^O linewidths. 
in th e range ~ 0.14 — 0.61 km s~^ (jMuench et al 



20071 ) ■ are typically 1.4 times broader than NH3 
which, like N2H+, traces only the gas ass ociated 
with t he dense compact core. Furthermore, [Onishi 
( 19991 ) found that the ^'^CO emission lines from 
the Pipe cloud were characterized by linewidths 
of ~ 1 km s~^. This leads in this case to a typical 
i^CO to NH3 linewidth ratio of ~ 2.8. 



The difference between the two scale factors 
may be caused, besides to the statistical signif- 
icance of the two samples of sources, also by 
real dynamical differences bet ween t he V ela-D 
and Pipe re gions. In fact, both lElial (|2007l ) and 
Massi et alj (J2007) find evidence for turbulent ac- 
tivity in Vela-D. In addition, the presence of mul- 
tiple IRAS embedded clusters, stellar jets and out- 
flows also suggest that overall Vela-D is in a more 
advanced evolutionary state compared to the Pipe 
nebula. 

We also not e that if we use the NH3 survey 
of IJiiina et all |l999l) and the "CO survey of 
Vilas-Boas et al.l (l2000V the average "CO to NH3 
linewidth ratios toward, e.g., the Ophiucus and 
Taurus clouds would be even larger than that es- 
timated for the Pipe cores. Therefore, the reader 
should consider the results obtained by applying 
the linewidth scale factors with caution. In the 
following, we will use the 1.6 scale factor obtained 
toward Vela-D as our baseline value, but we will 
occasionally discuss how the results would change 
if a larger (typical of other clouds) scale factor 
were used. 

2.3.3. Comparison with Previous Results 

Comparing Figures |3] and |4] we note that the 
core mass, Mcore — 3 — 8Mq, at which the crit- 
ical threshold Mcorc = -Mbe is crossed, is quite 
lower than the mass at which the cores appear 
to become gravitationally bound, which from Fig- 
ure |4] can be roughly estimated to be > 30 M©, 
when the scale factor for the linewidt hs is applied. 
Looking at Figure 12 of lOlmil (|2009l ). it would be 
tempting to infer that the characteristic mass at 
which the core mass function appears to break or 
depart from a power-law is similar with the criti- 
cal Bonnor-Ebert mass of Figure [3l However, we 
cannot compare these mass ranges because of the 
completeness limit. Thus, we cannot conf i rm for 
Vela-D the results obtained by iLada et al.l (|2008l ) 
for the Pipe cores, where the authors find a con- 
sistency between all of these critical mass ranges. 



As far as the results of lElial (|2007l ) are con- 
cerned, we have already seen that we associate 
a BLAST core to a single "C0(2 - 1) veloc- 
ity compon ent, whereas in the clump catalog of 
Elial (|2007l ) multiple gas clumps may be associ- 



ated with a group of dust cores, as detected by 
Massi et al.l ( 20071 ). Because of this and also be- 



cause Elial pOOTI ) do not attempt to correct the 
^^C0(2 — 1) linewidths, it is not surprising that 
their virial masses turn out to be a factor ^3 — 40 
larger than the estimated gas mass of each clump. 

On the other hand, assuming that either N2II+ 
or NII3 trace the bulk of the core emission, one 
can apply the 1.6 scale factor estimated in Sec- 
tion E12] to the "C0(2 - 1) linewidths. In this 
case we find that ~ 35 % of the Vela-D cores de- 
tected by BLAST wi thin t he ar ea covered by the 
"C0(2 - 1) map of JEhal (|2007t ) are near or over 
the threshold to become gravitationally bound, 
as shown in Figure |4l If we had used a larger 
linewidth scale factor, such as the one derived from 
the Pipe cores, nearly half of the Vela-D cores 
would result gravitationally bound, according to 
this analysis. However, the exact fraction of bound 
cores will be known only when the linewidth of 
each core will be individually measured with a 
high-density tracer. 

3. EFFECTS OF PRESSURE 
3.1. Cores Edge Pressure 

To further investigate the dynamical state of 
the Vela-D cores we now analyze the effects of 
thermal and non-thermal pressures on the cores. 
First, we compare the core mass with the esti- 
mated radius of the Bonnor-Ebert sphere, for a 
given edge pressure. The Bonnor-Ebert critical 
mass can be written in terms of its radius and of 
the pressure, Po , at the edge of the sphere (equal 
to the ambient pressure) as: 

MBE[AfQ] = 0.15 (Z^beM)' (Po/fc[cm-3K])i/2. 

(8) 
In the top panel of Figure [5] we have thus plot- 
ted for each core where a measurement of A/vir 
was available, the edge pressure Po/k, correspond- 
ing to a Bonnor-Ebert sphere (with given radius 
and temperature) in hydrostatic equilibrium, as a 
function of the Mcorc/AfeE ratio. In addition, we 
have also plotted the sources where Mcorc /A/vir > 
0.5. 

Given the small number of data points, the in- 
terpretation of the top panel of Figure [5] is dif- 
ficult. The number of data points becomes even 
smaller if we take into account the fact that the 
definition of Po is only valid for sources at or near 
the Bonnor-Ebert equilibrium. For this reason we 
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Fig. 5. — Top. Pressure, Po/k, at the edge of the 
Bonnor-Ebert sphere (equal to the ambient pres- 
sure) vs. Mcorc/M-QE, for aU points where a virial 
mass could be determined. The green symbols rep- 
resent cores with Mcorc > 0.5 Afvir (1-6 scale factor 
applied, see text; black symbols are for cores with 
Mcoie < 0.5Mvir). Bottom. Edge pressure calcu- 
lated assuming A/be = A/core (see text). Green 
and black symbols are as before. The dashed hor- 
izontal line corresponds to the ambient pressure in 
Vela-D, Pext/k '-5x10^ cm"^ K (see Section[ 



have drawn in Figure [5] two dashed vertical lines 
corresponding to Mcotc/Mbe = 0.5, 2 which en- 
close all sources having a mass within a factor of 
2 from the corresponding Bonnor-Ebert mass re- 
quired to be in hydrostatic equilibrium. For these 
points there is a wide distribution of Po/k values, 
from Po/k - 10^ cm^^ K to - 10^ cm^^ ^ how- 
ever, we tentatively note that most sources with 
Mcotb/Mbe < 1 have a higher (> 4 x 10'' cm"^ K) 
value of the edge pressure, compared to those ob- 
jects with Mcoib/Mbe > 1 (three sources only, 
with pressure < 3 x 10^ cm~'^ K). This is indeed 
what one might expect in cores with low density- 
contrast which are mainly confined by the exter- 
nal pressure, and not by self-gravity. Clearly, a 
larger source sample is needed to confirm that this 
trend is real. We also note that most sources with 
Mcoie/Mvir > 0.5 (with a 1.6 scale factor) also 
have Mcorc /-A/be > 1, a trend which is confirmed 
when a larger linewidth scale factor is used. 

In the bottom panel of Figure[5]we have instead 
plotted the edge pressure, Po/k, assuming Mbe = 
Afcoroj as a function of the A/corc/A/BE ratio, es- 
timated as described in Section [2?2] From Equa- 



tions (HI and ^ one can see that Po/k oc T^ Afg^. 
Therefore, for cores with Afcore/AfBE < 1, the 
edge pressure calculated as Po/k ex T'^ M^^^ rep- 
resents the external pressure required to make the 
core gravitationally unstable. On the other hand, 
in the cores where Afcoro/A/BE > 1 the quantity 
Po/k ex T'^ M~^^ represents the internal pressure 
required to make the core stable against gravita- 
tional collapse. 

Since the ambient pressure in the densest parts 
of the Vela-D region has been estimated to be 
Pext/k ~ 5 X 10^ cm-3 K (see SectionE^]) the bot- 
tom panel of Figure [5] is suggesting that the cores 
with A/corc/A/BE < 1 are actually likely to undergo 
gravitational collapse at some point in the future, 
unless an additional, non-thermal source of inter- 
nal support exist. This internal support could in 
fact be of turbulent origin, since most of the points 
with A/coic/A'/be < 1 also have Afcorc/Afvir < 0.5 
(see Section [2?3| . 

On the other hand, the cores with Afcoro/AfBE > 
1 show that the required internal, thermal pres- 
sure to remain in equilibrium is much less than the 
ambient pressure in Vela-D. Thus, these cores are 
also likely to undergo gravitational collapse, ex- 
cept possibly the objects with Afcoro/Afvir < 0.5, 



which have an additional source of internal turbu- 
lent support. 



3.2. Thermal and non-Thermal Pressure 

We now specifically consider the ratio of ther- 
mal to non-thermal pressure, Rp = (a/cTNT)^, 
for each core; here a = [/cT/(/x amu)] is the one- 
dimensional isothermal sound speed in the gas, 
ctnt — A^T/(81og2)^/^ represents the non- 
thermal component of the velocity dispersion, and 
AVnt - [AVi^aco " 81og2fcr/mi3co]'/', where 
AVi3co is the observed linewidth. 

In the top panel of Figure |6] we plot Rp vs. 
core mass with and without the scale factor for the 
linewidths fSection l2.3p . This figure shows that all 
of the cores in Vela-D for which we have a mea- 
sured ^'^C0(2 — 1) linewidth are below the Rp = 1 
line, and thus the non-thermal pressure would ap- 
pear to exceed the thermal pressure. Figure|6]also 
shows that the situation only marginally change 
if the linewidths are scaled with the 1.6 factor. 
Many more cores would have i?p > 1 if the larger 
2.8 s cale factor from th e Pipe cores were used. In 
fact, iLada et al.l (J2008I ) found in the Pipe cores 
that where both NII3 and C^^O are observed, the 
non-thermal pressure derived from the C^^O data 
is on average a factor of 4 larger than that mea- 
sured by the ammonia line. 

We then evaluate the average total internal 
pressure of each core in Vcla-D for which we have 
measured the ^"^00(2 — 1) linewidth. We note that 
this pressure will overestimate the pressure at the 
edge of the cores. Including both thermal and non- 
thermal contributions, the pressure is given by: 



^tot = P(a^ + crNT) 



(9) 



where p is the mean density of a core. In the 
bottom panel of Figure [HI we plot the quantity 
Ptot/k vs. core mass, where we note the high 
values of the pressure, due to the relatively high 
density of the cores and to the large ^'^C0(2 — 1) 
linewidths. These values are consist ent with the 
cores internal pressure estimated by iMassi et al.l 
(J2007[ ). who also found that the external pressure 
to the cores is of order Pext/k ~ 5 x 10^ cm""^ K. 
However, if the 1.6 scale factor for the linewidths 
is applied, the total internal pressure would vary, 
for most cores (Figure [HI), in the range Ptot/k ~ 
5 X 10^ — 4 X 10^ cm~'^ K, comparable to the am- 
bient pressures characterizing other SF regions. 
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Fig. 6. — Top. Ratio of thermal to non-thermal 
pressure for the Vela-D cores vs. core mass. Bot- 
tom. Total core internal gas pressure plotted as 
a function of core mass. Symbols are as in Fig- 
ure m The small green "-(-" signs correspond to 
the pressure calculated using a Larson-type scal- 
ing law (see text). 



In the lower panel of Figure [6] it is also evident 
the trend of increasing Ptot/k for more massive 
cores. We find a trend of increasing density with 
increasing core mass in Vela-D, but this trend is 
not enough to explain the increase in Ptot/k with 
core mass over more than two orders of magni- 
tude, as shown in Figured Therefore, the higher 
value of Ptot/k in more massive cores must also be 
a consequence of the increased thermal and non- 
thermal components of the internal pressure. In 
particular, the top panel of Figure |6] shows that 
the contribution of the non-thermal pressure be- 
comes dominant for high-mass cores. 

This is in fact consistent with the trend ob- 
served betw een pressure and mass when a simple 
scaling law (JLarsonlll98ll) is used between the line 
wi dth and th e densit y, AV ^ oc p~^/'^. According 



clusters in the V ela-D region (see lElial 120071 and 
Massi et al.ll201Clr ) . This would suggest that either 



to lAdams fc Fatuzzo (i 19961 ) this implies a scaling 
relation between the linewidth and the mass of the 
core, which can be written as: 



McoroiMo] - 466.5 (aoff[kms- 



(10) 



this relation can also be written using Eq. ([9]) for 
the "effective" sound speed, a^g = -Ptot/p, as: 



Ptot[dynecm-2] = 4.6 lOVfecm-^] (Afcorc[Me 



)l/2 
(11) 

which is plotted in the lower panel of Figure [6l in 
the usual units of Ptot/k. We note that the model 
pressure reproduces the trend followed by the data 
points, but for each specific point the pressure 
given by Eq. (jlip generally underestimates the 
pressure given by Eq. (jH]), though to a lesser extent 
when the scale factor for the linewidths is applied. 

3.3. What Supports the Vela-D Cores? 

We thus have two possible scenarios regarding 
the dynamical state of the Vela-D cores, if we con- 
sider the sub-sample of cores within the area cov- 
ered by the ^■^C0(2 — 1) map as a representative 
sample of the whole Vela-D. If we do not apply 
the scale factor for the linewidths, then almost all 
of the cores would appear to be gravitationally 
unbound (Figure |4]) and would have an internal 
pressure much higher than the core external pres- 
sure (Figure [S]). We note, however, that the aver- 
age value of the core deconvolved radius to sound 
speed ratio is (i?doc/a) = (3.2 ± 1.1) x 10^ yr, thus 
quite lower compared to the estimated dynamical 
age, ^ 1 to 2 X 10^ yr, of the filaments and stellar 



the gravitationally unbound cores are really tran- 
sient structures (e.g., p roduced by th e turbulence 
of the Vela-D region; see lEhal2007l and lMassi et all 
20071 ) and are thus constantly replenished, or they 



are externally supported. 

As for the latter case. Figure |6] shows that many 
of these cores cannot be pressure-confined. Then 
the question arises of how the gravitationally un- 
bound cores can be effectively confined against 
their own internal pressure to still form coher- 
ent structures. Because in this case self-gravity 
and surface pressure alone appear to be insuffi- 
cient to bind the cores, a plausible source of exter- 
nal confinement could come from the presence of 
static, helical magnetic fields within the filaments 
or clumps where the d ense cores have formed (see 
Fiege fc PudTit3l2000l ). 



We can estimate the magnitude of the field 
strengths needed to produce the external magnetic 
pressure which would add to the ambient gas pres- 
sure to compensate the cores internal pressure, i.e. 
Pb = B^ /{Sir) = Ptot — Pext- In this way we find 
that the required B would vary from a few tens 
to several hundreds of fiG. While these values are 
in agreement with measure d magnitud es of B in 
molecular clouds (see, e.g., lBasull2004l and refer- 
ences therein) , we would have to assume that all of 
this field is toroidal and constitute s the dominant 
component (JFiege fc Pudritj[2000l ). 

On the other hand, we know that the ^■^00(2 — 
1) linewidths are indeed overestimating the core 
internal velocity dispersion, as discussed in the 
previous sections. Thus, if we do apply the ap- 
proximate scale factor discussed in Section 12.31 
then ^ 30 — 50% (depending on the scale factor 
assumed) of the cores (mostly proto-stellar) would 
become gravitationally bound. The internal pres- 
sure of these cores still appears to be significantly 
higher than the pressure of the ambient medium, 
but because the cores would be confined by self- 
gravity the high surface pressures are of no concern 
for the overall long-term stability of the cores. 

The cores that are not gravitationally bound, 
and with internal pressures higher than the sur- 
face pressure, would still need a source of external 
confinement, unless they are transient structures, 
as discussed earlier. However, we find that many 
(~ 40%) of these unbound cores (both starless 
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and proto-stellar) would have an internal pressure 
Ptot/k ^ 2 X 10^, which is a factor of 4 higher than 
the estimated average external pressure in Vela-D, 
Pe^t/k - 5 X 10^ cni-3 K (JMassi et al.ll2007t) . thus 
partially reducing the imbalance between internal 
and external pressures. 

4. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the dynami- 
cal state of the population of the dense and cold 
cores detected by BLAST in the Vela-D molecu- 
lar cloud, utilizing the sensitive maps at 250, 350 
and 500 /im obtained by BLAST during its 2006 
LDB flight from Antarctica. We can summarize 
the results of this paper as follows, although some 
of our conclusions are still dependent on an exact 
measurement of the internal velocity dispersion of 
the cores: 

1. The median Mcotb/Mbe = 0.9 for all sources 
in Vela-D, i.e. very near the critical value, assum- 
ing thermal support only. If we restrict ourselves 
to mass ranges for which we are progressively more 
complete, our sources exceed their Bonnor-Ebert 
mass by a median factor of Mcore/-^BE = 1.5 — 3.3. 

2. If the BLAST cores were purely isother- 
mal structures, then almost all would be either 
confined by the external gas pressure, or by their 
gravity. According to the virial analysis of the 
cores with an unambiguous line-of-sight ^^C0(2- 
1) spectrum, almost none of them would appear 
to be gravitationally bound and would also be 
dominated by non-thermal pressure. However, we 
think this is likely a consequence of the "'^^00(2-1) 
linewidths overestimating the core internal veloc- 
ity dispersion. Despite the uncertainties involved 
in this procedure, we have attempted to correct for 
this effect and found that ~ 30 — 50% of the cores 
would then turn out to be gravitationally bound. 

3. The average total internal pressure of the 
cores appear to be higher than the ambient ex- 
ternal pressure. If most cores were indeed grav- 
itationally unbound then self-gravity and surface 
pressure would be insufficient to bind the cores 
and form coherent structures. After application of 
the scale factor to the ^■^C0(2-l) linewidths, a sig- 
nificant fraction of the cores (mostly proto-stellar) 
would actually become gravitationally bound, and 
thus they would not need to be pressure-confined. 

4. The cores that result gravitationally un- 



bound even after application of the scale factor 
to the ^'^C0(2-l) linewidths, would still need an 
additional source of external confinement. If the 
confinement mechanism were an external magnetic 
pressure, the required B would vary from a few 
tens to several hundreds of fj,G. Measurements of 
the BLAST cores velocity dispersion is necessary 
to allow a more accurate dynamical analysis. 

5. Alternatively, the gravitationally unbound 
cores could be transient structures, continuously 
produced in the region by turbulence. 
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